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AKSfl4ACT 


This  p.ipor  duiiib  with  the  itcr.uivu  solnt  Ion  ol  mm- 1 i near  eipi.it  ions  1 (x)  “ 0. 
integral  inlomation  on  i wtiich  is  p,ivcn  l>y  : ( ) , 1 ' f , . . . , 1 l<t)ilt. 

polaiory-integial  metliod  wliitli  uses  integrii  Inlormation  and  wliicli  has  maximal  order 
equal  to  s<-3.  Since  the  maximal  order  ol  iterations  wliicli  use  1 (x^)  , . . . , I Nx^j)  is 
additional  inlormation  given  hy  tiie  Integral  Kl)dt  increases  tlie  order  by  two. 


We  consider 
We  deline  an  inter 
ol  convergence 
equal  to  s+1,  tlie 


V 


0 


1 . ISTKODUmON 


We  consider  tite  solution  of  tlie  nonlinear  scalar  equation 


(1.1)  fix)  - 0, 

wliere  f Is  a complex  function  of  complex  variable. 

In  most  papers  which  deal  with  stationary  Iterative  methods  for  (1,1)  It  Is  assumed  we  know  the 
standard  inforniation  for  t (Wozniakowski 

wtiere  s «-■  I and  x^^  is  an  approximation  to  the  solution  cr.  Hie  maximal  order  of  convergence  of  such 
methods  is  equal  to  s + I (Traub  [6AJ,  Wozniakowski  [73)).  We  raise  the  questi'n  low  other  types  of 
inlormation  can  be  used  in  iterative  processes  and  what  is  *'he  maximal  order  ol  convergence  ior  this  in- 
format ion. 

This  paper  deals  wiiti  Integral  inlormation  which  consists  of  the  standard  inlormation  Ul^  and  addi- 
tionally the  value  ol  an  nlegral.  Tims 

’'o 

(1.2i  31  , - n(x„) i'®’(x.),  ■■  Id)  dt’ 

- 1 ,s  n ’ w . 

wtioru  Yq  is  a ctimplfx  mimht*r  defined  in  Section  J, 

in  Section  2 we  deline  an  1 nterpoi.Uory  - integral  method  1 , which  uses  integral  inlormation 
31  to  estimate  n and  in  Section  7 we  prove  its  order  lor  s • 1 is  maximal.  Sections  , 5 and  6 con- 

tain  theorems  aiiont  the  convergence  ol  l_j 

Wozniakowski  7A  j deiined  lor  tiie  generalized  inlonnution  '3  an  order  ol  infumiation  pCH)  and  proved 

it  is  equal  to  tlie  maxinal  order  ol  convergence.  in  Section  7 we  prove  lliat  for  s ?:  I and  lor  suitable 

^0 

cliosen  y , pC'l  ) - s + i.  Since  p('3  ) = s t I,  the  additional  ini  orra.it  ion  given  hy  fitldt  increas- 
('  -1,S  S 

es  the  or' er  ol  Inlormation  hv  two.  For  systems  ol  nonlinear  equations  similar  resiiils  can  lie  proved 
and  will  he  reported  to  a liitiire  paper. 


1 


2.  INTbKi'OiJVrOKY  - iNTC  jU.  IIKKAIIVT-;  MICTIlOi)  f_,  ^ 

Let  us  consider  the  solution  oi  the  nonlinear  equation. 


(2.1)  t (x;  - 0, 


wliere  f:D  - (T  , U Is  an  open  subset  of  (J]  , (£  denotes  the  set  of  complex  numbers.  Let  of  € D be  a 
almple  zero  ol  1,  Ho)  “ 0 / f'(o).  An  Interpolatory  - integral  method  is  defined  as  follows. 

Let  Xj  be  an  approximation  of  We  assume  that  the  Information  on  f is  given  by 

’'i 

(2.2)  ')t  , - , (x  ;f)  - ff(x  ) f^®Nx  ),  ^ f(t)  dtl, 

-I,  Si  1 1 y. 

where  y^  depends  on  x.,  f^'^^x^),  k - 0,1 y^  ^ x^  and  is  defined  in  Section  3.  If  s - 0 then  y^ 

can  depend  on  x^,  f(x_.),  x^  ,,  f(x^_,).  The  value  of  y^  will  be  chosen  to  maximize  the  order  of  itera- 
tion. 

The  information  consists  of  the  standard  information  given  by  f(Xj) f^*  (x^)  and  additionally 

the  value  of  the  integral.  Next,  let  w^  be  an  interpolatory  polynomial  of  degree  at  most  s +1  such 
that  ; 


(2.3)  W|'‘’(X()  - ^^''^x^) 


k “ 0,1,...,s, 


(2.‘»)  r Wj(t)dt“’  f(t)dt. 

■ ^i  >'i 

if  w^  exists  then  the  next  approximation  x^^j  in  1 j ^ method  is  defined  as  a zero  of  polynomial  w^^, 
(2.5)  + ‘ 0. 

with  a criterion  to  make  x.^^  unique.  We  shall  now  prove  that  w^  exists  and  is  unique.  Let 

X 

(2.0)  F(x)  = ' 1 (t ) dt. 

^1 

and  let  g.  be  a pt)lynomial  oi  degree  s -♦  2 such  that 
(2.7)  gj(y^)  - Ify.)  - 0, 

(2.0)  gj''’(x.)  = i'^'‘’(Xj),  k - 0,1 s+). 

Ttius,  g.  is  a llermite  interpolatory  polynomial  for  K.  The  assumption  y,  / Xj^  Implies  the  existence  and 
tile  uniqueness  oi  g.  Set 


(2.9)  w.  (x)  - g'.  (x  I . 
1 I 


2 


Tl\en  (2.3)  and  (2.4)  are  satislied  wliitli  coiiiplc-tes  the  pruol  . Moreover,  1 rom  (2.8)  and  (2.9)  follows  an 


error  formula, 


s+2 


(2.101  K(x'  - Rj(x)  - (x  - y.)(x  - X|)  ‘ i:j(x) 


where 

(2.11)  G (X)  - 0 (x,f)  - j‘  r’...  +t,(Xj-y.)  4 t^^3, 

‘ 0 0 0 


r 


(x  - x^))  dtj  ...  dt^^3* 


Dlf*'eren'latlng  (2.10)  we  get  from  (2.11), 


(2.12)  f(x)  - w^(x)  - R(x) 


where 

R(x)  R(x,f)  “ (x  - Xj)  {[  (s+2 ) (x  - y^^)  + x - x^]  C^(x)  + C.(x)(x  - x^)(x  - y^)!. 

3.  DKKINITION  OF  A LOWER  LIMIT  OF  THE  INIT.GKAI. 

We  want  to  define  yj^  to  maximize  the  order  of  l_j  Setting  x • or  In  (2.12)  we  have 

(3.1)  - w^(o)  “ R(a) 

Let  us  assume  for  a moment  that  w^  has  a zero  x.^j  sufficiently  close  to  a simple  zero  or. 

w (a,) 

X . 5 ■ . I + 0((x_,  - nr)')  “ 0(R(a)) 
i+1  1 + 1 

We  see  that  the  order  of  Iteration  depends  mainly  on  R(a).  Hierefore  we  shall  choose  y^  to  minimize 
R(a)  In  a certain  sense.  From  (2.12) 

(3.2)  R(a)  “ (a  - x ^ ) (s+2 ) (o  - y^)  + m - x,  J G,(a)  + GjforHo  - x^)(ot  - y^). 

As  Gj(a)  and  C|(a)  are  In  general  unknown  we  want  to  minimize 

(3.3)  max(|(s  + 2Ha  - y^)  + o - x^,,  | (a  - x.l  (a  - y,)|)- 

One  can  verify  that  the  minimal  value  ot  (3.5)  is  lor  yj  equal  to  y 


(3.4)  y = a + 


• X . 

1 

2 4 la  - x,| 


) 


As  we  do  not  know  a we  have  to  replace  It  by  an  approximation  to  a,  2^  which  depends  only  on  the  stan- 
dard infonsatlon.  - 2^(x^,f(x^) and  x^.  ,f  s - 0 then  z^  - (x^.,  . f (x^. , ) .x^ . f (x  ^)X 


We  define  y^  as 


(3.5) 


s + 2 + |zj  - x^l 


It  can  be  proved  that  one  can  drop  - xj  In  the  denoa.lnator  without  the  change  of  the  order. 
Finally,  y^  Is  defined  by 


(3.6)  y - z + 


^ - "l 


1 1 s + 2 • 


Hence,  from  (3.6)  and  (2,12)  we  get 


(3.7)  f(x)  - Wj(x)  • R(x) 


for 


R(x)  - K(x,f)  - (x  - x^)®"^’ ((s+3)(x  - z^)  G^(x)  +Cj(x)(x  - x^)  ^ ^ 

where  Is  the  Interpolatory  polynomial  defined  In  Section  2. 

A.  THE  CONVERGENCE  OF  THE  ITERATU'E  METHOD  I FOR  s 5?  1 

In  the  previous  section  we  have  seen  that  the  order  of  Iteration  mainly  depends  on  R(a).  From 


(3.7) 


D/  V / .s+1,  (s+3)(a-z  ) + X -a 

K(of)  - (a  - Xj)  f(s+3)(a  - Zj)  Cj(a)  + G|(a)(a  - Xj) ^ ^ 


s + 2 


Hence,  to  assure  the  maximal  order  of  lor  s > 1 It  suffices  to  define  approximation  z^  using 

Newton  method 


(4.1)  z^  - 2^(x^,f(x^) ■ x^  - ^ 


I (x,  ) 


(Xj) 


Ttieorem  1 

If  s 5;  1 and 


1 • 

• IS  a continuous  function  on  K(r,,K)  where 


I ■ 0,1,... 


f(a)  • 0 4 f’(a) 
Kla.R)  ■ Ix  : lx  - a| 


■-  Kl,  R = K(n  - max^tl)_cr  ^ ^ 


where  the  constant  C is  defined  below. 


4 


2.  a real  number  |~  0 Is  such  that 


2 M 

P"  li([~)  I and  — — f~  1 where 


'*1  " |f  ^\x)|  i - s+2,8+3; 


x€K(o,R) 


M,  ■ sup  |f(x)|;  C - C(n  - 2-^ fT 2 : 

‘ x€J  . 12  ,-)^ 


3.  Xq  € J where  J « [x  : |x  - cr|  - , 

then  the  sequen.e  [x^)  generated  In  I j ^ has  the  following  properties: 
(1)  X,  S J,  Vi  , 

c ^ 3 

(ll)  - a • ' 


|A  J < A Vi , and 

3 + 1 M a\(  , 

V,  I ' ' (s42>:  (s+4);  * s+2  ' 


(111)  iimx  -i>  and 

l-~- 


’‘i  + 1 ■ 

lim  — “ H where 

, , , s+  3 

i-~x  (x  - -y) 


It  - (-1) 


♦ 2 fr'(r>)  . l ^ 1 ^ 

(s+2):  f’(^Ms+2)j 


Proof  of  (1)  (liy  Induction). 

Let  us  assume  that  x^  f J.  From  (2.12)  w.fx)  ■ 0 Iff  x = lt(x)  where 


nr  ♦ ..  . K ( X ) 

Ii2i 


^ r7^ 


1 1 X / ot 


if  X = f» 


Now,  using  (■^.2),  n.7),  (2.11)  and  condition  2 one  can  verily  tti.  ► tire  assumptions  ol  the  Brouwer  fix- 
point  ttieorem  hold  tor  it  in  tl.e  set  J.  Hence,  tiiere  exists  x f J,  such  titat  Wj(x)  “ 0.  So  x.^j  f J 
and,  by  induction,  x.  • I i “ 0,1,... 

1‘rool  ol  ( i i 1 


from  (4.2)  and  (1.7)  we  get 


where  e^  “ Xj  - "i.  and  satisfies  tlie  relation 

2 

Zj  - O'  - C|(x^  - -r)  . 

An  upper  bound  on  A^  on<  can  llnd  using  an  assumption  2 and  (2.11). 

I'rool  ol  I 1 i i ) 

from  (■) . 1 1 and  ( I ) , 

|X(,,  - a|  •:  a|x,  - 

and  thus 

s+2  , 

i><i  + , - -vj  ' (A  P )‘+  jxji  - Vi  . 

from  assumption  2 lollows  that  A hence. 

Uni  X *=  fy, 

i-*a> 

finally  from  (4,5)  and  (4,4)  it  can  he  shown 


1 ini 


P+i  - ^ 


(x. 


s+3 


which  completes  the  proof  of  'riieorein  1. 


I im  A ^ 

J 


li, 


6 


In  general.  B Is  not  equal  zero  (see  point  (ill)  which  n«t.ns  that  s+3  is  the  order  of  the  Inter- 
poUtory  - integral  method  l.,^^  for  s • I.  (Trauh  [64).  Wo^niako«skl  [74]).  Note  that  iterative  meth- 
od. which  use  only  the  standard  Infonnatlon  l(x^) have  orders  at  most  sfl.  Additional  in- 

formation  given  by  ^ f(t)dt  Increases  the  order  of  1.,^^  by  iwo.  Tl,e  usage  of  the  I.,  method  in 
practice  is  profitable  if  the  evaluation  cost  of  the  value  of  Integral  is  approximately  equal  to  the 
evaluation  cost  of  function  or  its  derivatives. 

5.  THE  CONVERGENCE  OF  THE  f.^ 

Now,  we  assume  s - 0,  which  means  that  the  information  is  of  the  form: 


^ - "l 


-1.0 


(f(x^),  f(t)  dt)  where  Vj  " + 2 


Note  that  we  cannot  now  define  by 


the  Newton  method  as  we  do  not  know  the  value  of  the  first  deriva- 


tive. Let  be  given  now  by  the  secant  method. 


*1  “ '‘i  • f(xj  - f(x"“T  ^^"l^  • 


In  this  case,  the  interpoUtory  - Integral  method  I_,^„  is  a one-point  method  with  memory  (Traub  [64]), 


Theorem  2 


if 


1.  f^^^  is  a continuous  function  on  K(»,R)  - (x:  [x  - o|  < R),  where  f(y)  - 0 )<  f'(or) 

r) 


K - K( 


where  the  constant  C is  defined  below. 


2.  a real  number  ] 0 is  such  that 


2M, 


• h(r ) I and  ^ i 1, 


where 


for 


h(f)  = ^1(1  + cH  ■ + cf ) + 2rn 


inf  I ‘ I,  Mj  = s"P  If^'^^x)!  i = 2,3  ; 

' X - I ..  -L'/-.  U\ 


’ ' ‘ x.3<(',,R) 


m: 


1 


M;  = sup  |i"(x)l,  C - ( (|  ) 
xr.I 


7 


).  Xy,  Xj  ' .1  wlieru  .1  - (x;  |x  - a|  P ),  then  the  sequence  [x.},  generated  In  I ^ has  the  follow- 
ing properties; 

(1)  X;  ^ J 1 “ 0,1,... 

(11)  x.^,  - > - A,(Xj  - (x^_,  - 

whe  re  | A J • A 1 « 0 , 1 , . . . , 


, M,  M, 

A - ^ IdiCn  + C)  5^  + rf  [3(1  + cf)  + 2]1, 


( i i i)  1 im  X 

i-*CX. 


2 24 

and  moreover 


’‘i.i  - « 


lln  — 

l-“  (Xj  - i)^lXj_j 


— - B where  B - i f , 

o)  4U-((vW  » 


X - -E- 

ilvl  lira  ) — rjr^  ” where  p « 1 + ^ . | 

i->»  ‘’‘l  ■ 

The  proof  of  this  theorem  Is  omitted  since  it  Is  similar  to  the  proof  of  Theorem  1.  From  (iv)  follows 
that  1 is  the  order  of  the  intcrpolatory  - Integral  method  1 


.1,0- 


h.  IHJ.  CONl’BRGKNCK  OK  TIIK  I , METHOh  FOR  MirLTIPI.i;  ZEROS 

-I  ,s 

Let  us  now  assume  that  s 2 1 and  a is  an  m fold  zero,  l.e.. 


t(->)  • f'(a)  - ....  - f^'""'^a)  = 0 ^ f^"\a), 

where  m ■ s. 

Ttie  inlomatlon  is  given  hy 


21 


“ ff(x  ),  f(x  ) f^’’\x ,),  f(t)  dt) 

I ^ I 1 I » 


where 


z j - X, 

y.-  “ ' ~ ■ for  z » I (x,;f) 

1 U I s 1 


i i ,s  ^ 


The  notation  i^  ^ fx^;!)  Is  used  hy  Traiib  64  1 and  Wcznlakowskl  [73].  Tlie  character  of  convergence  of 
the  l_.  method  in  this  case  is  given  by  Theorem  3. 

“ I , s 


Theorem  3 

I ] s 1 and 


{ s + i ) 


••  1 is  a continuous  lunction  on  K(v,R)  = (x:  j x - aj  • K3  where 


8 


- R(P)  - ^ -+JT  ^ 


where  the  constant  D Is  defined  below. 


2.  a real  number  j~  0 is  such  that 


(Of'"  • ''<r ) 1. 


h(D-T^2  (r-Df")  (.*3)  (r  + Dr")-^2r  . 


- sup  |f<*-^’)(x)| 


V - Inf  , M - sup 

xfJ  (X  . o)""  x€J 


M ■ sup  |f  (x)  I 1 “ s+2,  s+3, 

‘ xrK(o.R) 


V (s+d: 

vm 


M , i S±i., 
V (8+d:  ^ 

m 


/M  , Y 2±L  £±i., 

...x  /^±j — \ 2 • r ’’ 

<“>  I 


3.  Xg  e J,  j “ {x:  |x  - ->1  ' P} 


(i)  Xj  f i • 0,1,... 


, llm  Xj  “ 
i -♦'x 


(iO  Ixj^,  - '■<1  '■  [xj  - >1 


S*^  1 

where  p “ min(-— ,2). 


Moreover 


Aj  • A, 


..1  ( s+1 

s-H  I -p  M 

m ^ 


2 (s4-3UHDr 

s+2 


9 


r 


(111)  Urn 


S+1+P 


s H where 


m/  ml 
/_\ 


(m) 


(a)  1 

\ 2±L  . 

If 

1 f 0 

If 

1 - 0 

• n,l\ 

— - P,0\lf^"’Na)|-(s+l):/ 


m 


(sf2); 


+ c 


2-p,0  (s+4) 1 s + 2 


^1,0 

(Iv)  p(m)  " (s+l+p)/m  Is  the  order  ot  convergence  of  I , method  for  m multiple  zeros. 

- 1 , s 

The  proof  of  this  theorem  Is  omitted  since  It  Is  similar  to  the  proof  of  Theorem  1. 

7.  HAXLMALm’  OF  I , 

-1  ,s 

Let  t , be  a class  of  stationary  Iterative  methods  tp  , which  use  Information  , and  which 
-1,s  -1,8  -1,s 

have  well  defined  order  p((p  ) (Wtznlakowskl  [7A]).  From  Theorem  1 It  follows  that  the  Interpolatory- 

-1  ,s 

Integral  method  I ^ ^ belongs  to  s 2 1. 

Now  we  shall  prove  that  I , has  maximal  order  In  the  class  Y , , l.e., 

p(I  . ) « sup  P('P  , ). 

’ ’-l.s^-l.s  ’ 

Woznlakowskl  [74]  defined  (Definition  7)  the  order  of  Information  p(!1’)  and  proved  It  Is  equal  to  the 
maximal  order.  Thus,  It  suffices  to  show  that  In  our  case 

p(I_,,^)  - P(7>.,^3).. 

Let  ^ f(t)  dtl,  for  any  y^^  « yj(x.,f(x^) f^®^(Xj^)). 


Theorem  4 


linen 


p(51_,  s + 3. 


z^  - X^  f(x.) 

If  s ^ 1 and  y^  ■ ^ ■ ’'i  - 


then  p(71  , ) “ s + 3. 

- I , s 


Proof 


Let  ^ be  the  class  of  complex  functions  of  complex  variable  which  have  a simple  zero  and  which 
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I 


are  analytic  In  the  neighborhood  of  o (see  Definition  I in  Woznlakowski  [74]).  We  recall  the  definition 
of  the  order  of  information.  Let  f € and  {t^lc  where 

i\a)  " 0 

(7.1)  f^(a.)  - 0,  1-0,1 11m  a - ft, 

l-«  ‘ 


(7.2) 


lint 

i-KO 


for  k - 0.1 , 


>8 


cS. 


g(a)  " 0. 


Next  let  us  assume  that 


I 


t 


► 


(7.3)  j(x^,f) 


''-l.s<’‘l-'l> 


Vi  , 


Where  (xj  is  an  arbitrary  sequence  converging  to  a.  Ut  be  an  Interpolatory  polyno«lal  of  degree  at 
most  s + 1 defined  as  follows: 

’’-).s(’‘l*“l)  - Vi. 

Thus 

f(x)  - f^(x)  • f(x)  - w^(x)  + w^(x)  - fj(x)  - R(x,f)  . R(s,fp. 

From  It  and  from  (2.12)  It  follows 

(7.4)  ft  - ar^  - 0((f(ft)  . f^(ft);)  - 0(|R(ft,f)|  + |R(ft,f^)|)  . 

- 0(|ft  - . |(s.2)(ft  - y^)  + ft  . xj  -f  (ft  - . (ft  - yj) 


Moreover,  we  shall  show  that  Mus  bound  Is  sharp.  Let  C be  a number  defim 

ft  - y 

I if  11m 

I-HB 

c 

Let 


I d as  follows : 


1 

s+2 


where 


h.(x)  “ (x  - Xj) 


s+l+C 


(X  - bj)  1 - 0,1,... 


Setting 


(s+2+C)  y^  + Xj 

s + 3 + C • 


(7.5)  f^(x)  - f(x)  + hj(x),  Vi 

for  any  function  f J , (K^,)  . o,  fj(a.)  “ 0)  one  can  verify  that  conditions  (7.1),  (7.2)  and  (7.3) 
hold.  Next,  there  exist  constants  C.  0 s.icli  that  1 ini  C - C 0 and 

1 i-t£  1 
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1 


i 


s+I+C 


I (s+2+t)  (o-yp+a-Xj^l 


(7.6)  |o  - >j|  - Cj|f(af)  - fj(ar)|  - Cj|or  - Xj 


I"  ■ M - i^sc 


O - X 


wiiich  proves  that  (7,4)  is  sharp. 

Krom  (7.4)  and  (7.6)  it  follows  Chat  lor  any  y the  order  of  infomation  p = p(!Jl  ) exists.  Let 

“ I , s 

US  assume  tliat  p s+3.  Let  e 0 be  a number  such  that  p-e  s+3+e.  For  f and  [f^  1 given  by 
(7.5)  we  get  from  (7.6): 


+ “ lim  sup 


|xj  - orl'*'*  i 


lim  sup 


|x^  - O' 


s+3-re 


which  is  a contradiction.  Hence  p ^ s+3  for  any  y^  - y^(x^ ,f (x^) . . . f (x^) . Now  we  shall  show  that 
the  above  estimation  of  p is  achievable  for  sit. 

Indeed,  setting 


* ’'l 

^1  “ ^ ^ ■ ’‘i 


we  have 


lim 


Of  - y 


l-<»  “ ■ ’'l 


i - J_  / r 

- r_T  <s°  C 


It  is  easy  to  ^»rify  that  from  (7.4)  and  (7.6)  It  follows 


f(xp 

f'(Xj) 


1). 


P<^I,s'  ■ ® + 3* 

which  completes  tl.e  proof  of  Theorem  3.  ■ 

Krom  Theorems  1 and  3 we  get 

Corollary  1 

The  interpolatory  - integral  method  1 , is  maximal,  i.e.. 

- 1 , S * ' 
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